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We recently shown, using tight-binding calculations, that 
nonequilibrium valley polarization can be realized in 
graphene, when the current is injected through "val- 
ley filter": a ballistic point contact with zigzag edges. 
Here we demonstrate, that the effect is surprisingly ro- 
bust against changing the crystallographic orientation 
of the filter axis. Namely, the output current remains 
polarized unless a point contact has perfect armchair 
edges, at which two subblattices are equally represented. 
The polarization is inverted when the filter orientation 
crosses the amchair line and, subsequently, dominating 
subblattice index of terminal atoms changes. In a bended 
graphene strip, the valley-polarized current can be di- 
rected towards one edge providing a theoretical possi- 
bility to observe a zero-magnetic-field analogue of the 
well-known Hall effect. For the valley valve, build of two 
valley filters in series and controlled elecrostatically by 
a gate voltage, the conductance-to-gate characteristic is 
inverted when 7r/3 vertex is placed between two filters. 
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Generic valley filter attached to heavily-doped rectangular 
graphene leads (shadow areas). Notice the constriction edges 
are neither perfect zigzag nor armchair. 
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1 Introduction Graphene, the new and probably the 
most intriguing one of carbon allotropes, have become a 
subject of intense research focus in last few years. As its 
low-energy spectrum is described by the Dirac-Weyl Hamil- 
tonian of massless spin-1/2 fermions graphene of- 
fers the unique possibility to test predictions of relativis- 
tic quantum mechanics in a condensed phase. However, in 
graphene, Dirac electrons have an additional degeneracy, 
corresponding to the presence of two different valleys in 
the band structure. This phenomenon, known as "fermion 
doubling", makes difficult to probe the intrinsic physics of 
a single valley in experiments, because usually the contri- 
butions of the two valleys to a measurable quantity exactly 
cancels each other. An interesting example is production of 



a pseudo-magnetic field, which affects the electron in a sin- 
gle valley, by a lattice defect or distortion ||2l. The field has 
an opposite sign in the other valley, so the effect remains 
a playground for theoreticians as long as both valleys are 
equally populated in an experiment. These are the funda- 
mental reasons to look for a feasible and controlled way of 
the valley-isospin manipulation in graphene. 

The practical consequences of eliminating fermion dou- 
bling effects are even more exciting: The spin-based quan- 
tum computing in graphene quantum dots |3|, exploiting 
the superior spin coherence expected in carbon nanostruc- 
tures, is one of the most recent examples. More gener- 
ally, the potential of graphene for future electronics rests 
on the possibility to create devices that have no analogue 
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in silicon-based electronics H. Briefly speaking, it seems 
that valley degree of freedom might be used to control an 
electronic device, in much the same way as the electron 
spin is used in spintronics 131 or quantum computing 0. 

Therefore, a controllable way of occupying a single 
valley in graphene (in other words: producing a valley po- 
larization) would be a key ingredient for the so-called "val- 
leytronics". The author, Tworzydlo, and Beenakker recently 
proposed an electrostatic method of valley polarization in 
a single-mode quantum point contact with zigzag edges 
i7i. The work was build on earlier findings |8ll9l [T0l[m 
[T2l[T3 . 14 1 for graphene ribbons (long and narrow ballistic 
strips), that they may support a propagating mode arbitrar- 
ily close to the Dirac point, and that this mode lacks the 
valley degeneracy of modes that propagate at higher ener- 
gies. For armchair edges of the ribbon, this lowest propa- 
gating mode is constructed from states in both valleys, but 
for zigzag edges only a single valley contributes ll8ll9l[T0l 
fTTlfTZ l. In accord with time-reversal symmetry, the mode 
switches from one valley to the other upon changing the 
direction of propagation. We found |7 |, that depending on 
whether the Fermi level in the point contact lies in the con- 
duction or valence band, the transmitted electrons occupy 
states in one or the other valley of the band structure. More- 
over, the effect is not limited to long and narrow ribbons, 
ia fact the filter of the same length as width produces a po- 
larization better than 95%. It was also shown in Ref. Q 
that two adjacent valley filters function as a highly effec- 
tive valley valve, passing or blocking the current depending 
on whether the two filters have the same or opposite polar- 
ity. Valley-filtering point contacts may also affect magne- 
toconductance oscillations in a multimode graphene ring, 
showing a suppression of the lowest harmonic for opposite 
polarity ifBll . 

Here we complement our previous study Q by numer- 
ical simulation of the electron transport through a constric- 
tion with generic edges, evolving from zigzag to armchair 
upon "rotation" of crystallographic orientation of the sys- 
tem axis. The results show that the output current lies pre- 
dominantly in a single valley (chosen by a dominating frac- 
tion of zigzag edges), unless the edges are perfectly of arm- 
chair type. A basic illustration of the valley-polarized cur- 
rent behavior in a bended-graphene strip is also provided. 
Finally, we show the function of valley filter (and valve) 
may be inverted by an appropriate edges reorientation. 

2 Preliminary: Ballistic transport in graphene In 

this section, we briefly recall the basic theoretical and ex- 
perimental results concerning phase-coherent transport in 
graphene at the Dirac point. First, we review the theory of 
Ref. 1 14 1 by Tworzydlo et ai, inspired by an insightful pa- 
per of Katsnelson |fT6l , who used the Landauer transmis- 
sion formula to obtain the quantum-limited conductivity. 
Following the same approach, Tworzydlo et al. calculated 
the transmission probabilities of Dirac particles through 



a rectangular strip of W width and L length, attached to 
heavily-doped graphene leads. 

The low-energy excitations in graphene are govern by 
the massless Dirac equation for the four-component wave 
function ^ = {^a,^b,^'ai ^'b) *e general case (with 
Wx, 'I'x the amplitudes for different valleys, and X = 
j4, B the sublattice index). In the presence of external elec- 
trostatic potential V{y) the wave equation may be written 
as 



hv ( a^dx + cTydy 








If + eVir ^ e^, (1) 



where the velocity v = ^V^ra/h w 10^ m/s (with a = 
0.246 nm the lattice constant and r — 3.0 eV the nearest- 
neighbor hopping) plays a role of the speed of light, cr, are 
the Pauli 2x2 matrices, e is the electron charge, and e is the 
energy. We suppose the electrostatic potential V ^ Oin the 
strip {0 < X < L) and V = Voc elsewhere. In particular, 
the scattering problem for Eq. ([T]i can be solved analyti- 
cally for two different classes of boundary conditions iflTl : 

1 .Infinite mass boundary condition 
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with = for the metallic armchair, and = ±27r/3 
the semiconducting armchair (depending whether the 
strip width W is, or is not an integer multiple of three 
unit cells). 

The results for conductance in a heavily-doped lead 
limit 1^00! 00, including spin and valley degeneracies, 
may be summarized by 111 81 
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with a = 1/2 for infinite mass boundary conditions, a — 
for the metallic armchair edge, and a — 1/3 for the semi- 
conducting armchair edge. In the limit of wide and short 
strip W ^ L the formula (j6]) leads to the universal value 
of the conductivity 



= G 



(7) 



regardless the boundary conditions. The result (j?]) for the 
conductivity agrees with other theoretical calculations lfT2l 
[191 which start from an unbounded disordered system and 
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Figure 1 Theory and experiment for a ballistic graphene monolayer. Left panel: conductivity at the Dirac point as a 
function of the aspect ratio of a graphene strip. The curves are calculated by solving the Dirac equation for different 
boundary conditions: infinite mass boudary condition (black tick line), semiconducting armchair edge (red thin line), and 
metallic armchair edge (blue dashed line). The datapoints for the zigzag edge (triangles) are the result of a numerical 
solution of the tight-binding model ||23]| . Right panel: Experimental results of Miao et al. (red solid squares) compared 
with theory ll25l for metallic armchair edge (blue dashed hne). 



then take the limit of infinite mean free path I. For a long 
time, existing experiments [20,21J were quasi-ballistic, with 
I W ^ L, and finding ct k Ae^/h. The "missing I/tt 
factor" was found very recently |l22l, terminating the main 
controversy on ballistic graphene conductance. 

Ref. 1 22 1 provides also first experimental insight into 
the conductivity dependence on the aspect ratio W/ L. In 
Fig. [T] we compare the analytical results lfT4ll . reviewed in 
this section, with a computer simulation of electron trans- 
port for zigzag edge (where the analytical solution is lack- 
ing), and with the experimental results of Miao et al. l22l . 
First, we found (23 \ that formula ^ with a = (metal- 
lic armchair) match precisely the numerical data for zigzag 
edges (see left panel of Fig. [T] blue dashed line and tri- 
angles, respectively). This is caused by the fact that there 
is a nondegenerated propagating mode at zero energy in 
both cases. (Notice also the results for the semiconduct- 
ing armchair are close to that for infinite mass boundaries, 
since there is no propagating mode there.) This finding 
strengthen the theory significantly, as the presence of metal- 
lic armchair edge requires the width to be an integer mul- 
tiplicity of three unit cell and thus seems unrealistic in an 
experimental sample. Moreover, the recent findings 1241 
show that the zigzag boundary condition applies generi- 
cally for any angle 6* 7^ (mod tt/3) of the boundary 
(where 6 = labels the armchair orientation). 

For the purposes of this paper, conductivity dependence 
on aspect ratio in the range W/ L < 2 is particularly im- 
portant, as it indicates the degeneracy of the lowest prop- 
agating mode. The experimental results of Miao et al. still 
lies above the theory |25| in this range, with a relatively 
large disperion of datapoints (see right panel of Fig. [T}. 
The probable explanation is related to an uncertainity in 



W and L for small samples, which may cause a system- 
atic horizontal shift of the data. (Notice the conductivity 
increase rapidly when lowering W/ L below the threshold 
value W/L w 2.) One may also attribute the deviation be- 
tween theory and experiment to the smooth (long-range) 
disorder, which was shown |26 | to be responsible for en- 
hanced condictivity in previous experiments. Albeit suc- 
cesfully eliminated in work |22l for short and wide sam- 
ples, the smooth disorder may still play some role for long 
and narrow samples, because the charge screening become 
less effective in such quasi-one-dimensional systems. Nev- 
ertheless, the remarkable experimental progress achieved 
by Miao et al. allows us to believe, that ballistic graphene 
samples with W/L < 2 and electric transport dominated 
by a single, valley-polarized mode, are to be fabricated 
shortly. The aim of this paper is computer modelling of the 
measurements directly demonstrating the valley-isospin ma- 
nipulation in such nanostructures. 

3 Generic valley filter in graphene Here we con- 
sider the honeycomb lattice of carbon atoms in a strip con- 
taining a short and narrow constriction with quantized con- 
ductance. We suppose the constriction edges are parallel, 
namely, the narrow region is terminated at the top and the 
bottom along the same crystallographic axis nai + ma2 
(with ai, a2 the basis vectors). Therefore, the edges are de- 
fined by specifying a pair of integers (n, m) and the con- 
striction width W. The different ways of terminating the 
honeycomb lattice are illustrated in Fig. |2] (Fig. [3] shows 
two basic examples of constrictions with (1, 0) zigzag and 
(1,1) armchair edges.) For a long constriction, where the 
transport become one-dimensional, the crystallographic ori- 
entation of edges (n, m) determines the direction of prop- 
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zigzag edge (/i,m) = (1,0) 




Figure 2 Schematic illustration of the mechanism of 
valley filtering in graphene nanoconstrictions. Top panel: 
Honeycomb lattice of a carbon monolayer with ai = (a, 0) 

and a2 = (i, basis vectors (a is the lattice spac- 
ing). Generic edge of a finite system is defined by a pair 
of integers (n, m). Middle panel: Hexagonal first Brillouin 
zone of graphene. The conically shaped valleys touch six 
corners of the Brillouin zone, the three marked by black 
(white) dots are equivalent. Bottom panel: Dispersion rela- 
tion for the constriction with (1,0) zigzag edges. For the 
lowest mode, there is one-to-one correspondence between 
the direction of propagation (indicated by arrows) and the 
valley isospin. The polarity switches (i) when changes the 
sign of chemical potential /i, or (ii) when rotates the con- 
striction axis by tt/3 (the (0, 1) zigzag case). The valleys 
are equally mixed for armchair edges. 




X 



Figure 3 Schematic diagrams of valley filters with (1, 0) 
zigzag edges (a), (1,1) armchair edges (b), and the corre- 
sponding potential profile V{x) (c). Here, Woo = 17\/3a, 
W — 6\/3a, L — 12a. (Width of the constriction with 
armchair edges is W' = W cos 7r/6 = 9a.) Constrictions 
with predominantly zigzag edges allows one to control the 
output valley polarization by varying the electrostatic po- 
tential Uq. 

agation in the first Brillouin zone (see Fig.|2]i. As the low- 
est propagating mode in such a nanoribbon (i) lack the 
twofold degeneracy of higher modes, and {ii) may be ar- 
bitrary close to the Dirac point li8ll9l[T0l[m[T2l[T3l[T4l . the 
crystallographic orientation (n, m) also determines the val- 
ley polarization of current passing the constriction. 
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Figure 4 Operation of the valley filter with perfect zigzag edges (n, m) = (1, 0), fixed width W = 20\/3 a, and varying 
length L. Conductance (a) and valley polarization (b) for L ^ 36 a (solid blue lines) and i = 16 a (dashed red lines) 
as a function of /io/^ (with A — ttt/AO). The parameter /iq = /ioo ~ Uq is varied by varying Uq at fixed Hoo- (c) The 
polarization at fixed /iq = ±0.05 r w ±2Z\/3 (solid red and open blue symbols, respectively) as a function the filter length 
L. Notice the saturation of P for L > W. 



Our analysis starts from the tight-binding model of gra- 
phene, with Hamiltonian 



H = Y,n,\^){J\+J2^\i){^\ 



(8) 



The hopping matrix element r^j = — r if the orbitals \i) 
and \j) are nearest neighbors on the honeycomb lattice, 
otherwise Tij = 0. K: V{xi) varies only along the main 
axis of a strip (see Fig.|3]l. Namely, the potential equals Uq 
inside the constriction (x < < L, where L is the con- 
striction length) and zero everywhere else. By varying Uq 
at a fixed Fermi energy /ioo in the external leads (wide re- 
gions), we can vary the Fermi energy iiq = /ioo — Uq inside 
the constriction. 

The wide regions in Fig. [3] have (1,0) zigzag edges, 
and therefore support 2N + 1 propagating modes at the 



Fermi energy ^oo, which form a basis for the transmis- 
sion matrix t. Modes I — 1, 2, ... TV lie in the first val- 
ley [with longitudinal wave vector kj.a G (tt, 2tt)] while 
modes I = — 1, —2, ... — lie in the second valley [with 
kxQ € (0, tt)]. The zeroth mode / = lies in a single val- 
ley fixed by the direction of propagation. The conductance 
of the constriction is determined by the Landauer formula 



G 



2e2 



N 



N 



(9) 



-N 



The valley polarization of the transmitted current is quan- 
tified by 



P 



(10) 
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Figure 5 Valley polarization as a function of ^o/^ and crystallographic orientation (n, to) of the filter axis (specified for 
each curve on the plot). 



where we consider the case > 0, such that the zeroth 
mode lies in the first valley. The polarization P e [—1, 1], 
with P = 1 if the transmitted current lies fully in the first 
valley and P — —1 if it lies fully in the second valley. 
Throughout the paper, we took /ioo — t/3 to work with 
heavily doped graphene leads, while remaining at suffi- 
ciently small Fermi energy that the linearity of the disper- 
sion relation holds reasonably well. 

We have calculated the transmission matrix numeri- 
cally by adapting to the honeycomb lattice the method de- 
veloped by Ando for a square lattice ll27l . Results are shown 
in Figs. |4] and |5] We have fixed the width of the wide re- 
gions at Woo = 70-\/3 a, corresponding to 2N +1 = 29 
propagating modes. The narrow region has width W = 
20-\/3a (measured along y-axis, see Fig. [3^,b). We mea- 
sure the electrochemical potential /ip in the narrow region 
in units of the subband splitting for (1,0) zigzag edge strip 

A = \^^Ji'nTalW = -KhvjW. (11) 

For our parameters A = 7rT/40. 

In Fig. p] we demonstrate the operation of the valley 
filter with (1,0) zigzag edges. Figs.|4^ and|4|5 shows the 
conductance and valley polarization, both as a function of 
the electrochemical potential /ip in the narrow region. For 
positive /iQ the current flows entirely within the conduction 
band, and we obtain plateaus of quantized conductance at 
odd multiples of 2e^//i (as predicted by Peres et al. IIT2ll ). 
The plateaus in the conductance at G = (2n + 1) x 2e^ /h 
correspond to plateaus in the valley polarization at P = 
\/{2n +1). For negative /io the current makes a transi- 
tion from the conduction band in the wide regions to the 
valence band in the narrow region. The interband transi- 



tion destroys the conductance quantization — except on 
the first plateau, which remains quite flat in the entire in- 
terval — 3Z\/2 < /io < 3Z\/2. The resonances at negative 
/lo are due to quasi-bound states in the valence band |28] 
29 1 . The polarity of the valley filter is inverted for negative 
/iQ, almost without loosing the quality. 

The data for quadratic filter L = 36a « W (blue solid 
curves in Figs. and |^) attached to rectangular leads 
(cf. Fig. |3^) are almost identical with those obtained for 
the filter with triangular leads, studied in Ref. |7|, show- 
ing the lead-effects are negligible as the transport through 
the device is determined by a constriction. The new insight 
into the nature of valley polarization is provided with the 
data for shorter filters. For L ~ 16a w W/2 (red dashed 
curves) the polarization is still above 70% for |/io| ^ A 
(in comparizon to 95% for L w W), although the conduc- 
tance plateaus are smeared out. This is because the con- 
tribution to the conductance from evanescent and propa- 
gating modes become coparable for W/ L > 2 (notice the 
conductivity in Fig.[T]is not far from the universal limit in 
this range). In contrast, the contributions from evenescent 
modes in two valleys cancel each other in Eq. ( fTO] ), ex- 
paining the robustness of this quantity when changing the 
aspect ratio. (The fast saturation of the polarization with in- 
creasing L is illustrated in Fig.|4}; for two selected values of 
the chemical potential /io = ±0.05 r « ±.2A/3.) There- 
fore, it would be very benefitial to find a quantity, which is 
proportional to P and, simultaneously, feasible to measure 
in an experiment. One of the candidates is the transverse 
Hall voltage induced by a pseudo-magnetic field in bended 
graphene strip IS), briefly discussed in the next section. 
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Figure 6 Probability distribution |!Z'out(x)P for the electron after passing the valley filter (not shown) in a bended graphene 
strip. Left and right panel show the cases of an opposite chemical potential inside the constriction /io = ±0.05 r. The 
remaining filter parameters are W = 20\/3a, L = 34a, and {n, m ) — (1,0). The coordinate system of Fig.jsjis also 
shown. 



But first, we report the principal new results presented 
in Fig.|5] where we show the evolution of valley polariza- 
tion when changing the filter crystallographic orientation 
(specifies by (n,m) indices). For m/n < 1 the polariza- 
tion is almost unaffected for /io > 0, and changes gradually 
with increasing m/n for /ig < 0. Again, the polarization 
is more than 70% in wide parameter range. This is because 
the large Fermi wavelength at small /to causes the quality 
of the valley filter is quite robust against edge imperfec- 
tions, like a nonzero fraction of armchairs separting the 
perfect zigzags. For n — m (the perfect armchair edges) 
we obtain P w 0, as the valleys are equally almost mixed 
in this case. For (0, 1) zigzag edge case the polarization is 
inverted, as it corresponds to the rotation of crystalographic 
axis by vr/S. The filter behavior is therefore consistent with 
the recent theory of Ref. |24|, which has shown that the 
zigzag boundary condition applies generically to any crys- 
tallographic orientation of the lattice except the exact arm- 
chair orientation. 

4 Valley polarized current in bended graphene 

To illustrate the effect of pseudo-magnetic field |2| in a 
bended graphene strip we modify the off-diagonal part in 
Hamiltonian (|8]l. The vertical hopping elements varies now 
with y position across the strip (see Fig.jsj, namely, for \i) 
and \ j) the nearest neighbors 

-r(l + Oyi), if x.i = xj ^j^) 

This corresponds to the uniform perpendicular pseudo-mag- 
netic field (we use the same pseudo-field in both leads 
as well as inside the constriction). Here, we took OWoo — 
0.5, with = UQ^/ia. 



In Fig. |6]we plot the probability distribution in the out- 
put lead |tf/out(x) P (with x = {x, y) the position on a hon- 
eycomb lattice ||30J ), which is summarized over all incom- 
ming modes, each one contributing with 



f<!ul(x)|^ = 



y^4.;(/>fc(x) 



(13) 



kk' 



4^fc'i0fc(x)</'fe'(x), 



where (pki^) is the wave function of ^-th mode (note we 
consider an identical input and output leads). The compact 
formula for |if'out(x) P reads 



l'^out(x)|^ 



^«(x)P = Tr(ttMx)], 



(14) 



where p(x) is the density operator with matrix elements 

[p(x)]fcfc' = 0^(x)0fc'(x). 

The data presented in Fig. |6] shows the current leav- 
ing valley filter forms a well colimated beam, particulary 
in the case of positive polarization (left panel), when the 
electron stays in a conductions band when passing the con- 
striction. (The cyclotronic length in both cases is I0 = 
Moo/ {V^T0) « 100a.) The effect of bending is clearly op- 
posite for the two valleys, suggesting a nonzero Hall volt- 
age, proportional to the valley polarization P in the lowest 
order, may appear at zero field in a bended graphene strip 
containing the constriction. 

5 Valley valve with a vertex Following the idea of 
Ref. Q, we put now to different valley filters series, to 
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Figure 7 Different valley valves with zigzag edges (a-c) and their conductance (d-f) at fixed /ii = 0.05 r w 2Z\/3 as a 
function of [12! Green dashed lines at diagrams (a-c) mark the electrochemical potential steps between values indicated 
at the bottom. Blue solid and red dashed curves at conductance plots (d-f) correspond \a L ~ 12 a and 36 a, respectively. 
Green vertical line show the case [L\ = /i2 of no potential step between two valley filters building each valley valve. 
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show the conductance of the system is G w 2e^//i or w 0, 
depending whether the two filters have the same or oppo- 
site polarity. We analyze three systems, depicted schemat- 
ically in Fig. |7] first, build of two identical filters with 
zigzag edges (see Fig. |7^), second with crystallographic 
axis of one filter rotated by 7r/3 (see Fig. uh), and third, 
consists of two rotated filters building a nanoribbon with 
2tt/3 vertex in center (see Fig.jTj;). Thus, the current through 
each of these valley valves have to change its direction by 
9 = 0, 7r/3, and 27r/3 (respectively) when passing the 
center of a constriction. We use the lead width W^o = 
140%/3a (twice as large as in Section [s] to make sure the 
constriction is placed far away from lead edges) correspond- 
ing to 2N + 1 = 59 propagating modes in wide regions. 
The constriction width is Vt^ = 20\/3a (defined along 
the transverse direction, as in Fig. |3]l. Again, we perform 
the computations for two values of the constriction length 
L/a — 36 and 72. (So the valve is always build of two 
valley filters selected from those analyzed in Section[3]) 

The results are shown in Figs.|7|J-f, where we plot the 
conductance at fixed chemical potential in the first filter 
/ii = 0.05 r as a function of the chemical potential in the 
second filter ^2- The striking feature is an almost perfect 
inversion of the valve function each time when changing 
the relative angle 9 between crystallographic axes of the 
two filters in the steps of 7r/3 (with some loose of the qual- 
ity for /i2 < in = 27r/3 case, where the complicated 
romboidal shape of a hole-doped quantum dot produce a 
relatively dense spectrum of resonant states). In particular, 
for the situation without a potential step between the filters 
Ml = = 0.05 r (indicated by green solid line along the 
right panel of Fig. I?]) the current is passed for 9 0, 27r/3 
and blocked for tt/3. In other words, the valve function is 
reduced in this case to the mechanism of current blocking 
at a vertex of a graphene hexagon with zigzag edges, dis- 
cussed in Ref. ifTSll . 

6 Conclusions We have analyzed the operation on a 
valley filters with different crystallographic orientation on 
a honeycomb lattice. In a general case, the filter behavior is 
dominated by a zigzag line that is closest to its actual edge 
direction. The function of a filter is inverted when pass- 
ing the armchair orientation. The effect of pseudomagnetic 
field on a valley-polarized current in bended graphene have 
been also demonstarted. 

The conductance calculation for a valley valve shows 
it remains highly effective when changing the crystallo- 
graphic orientation of the one of two filters in series. The 
valve function is insensitive on whether the same (or oppo- 
site) polarity of the filters is achieved electrostatically or by 
changing the geometry. These findings constitutes an alter- 
native line of approach towards the practical realization of 
key ingredients for valley tronics. 
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